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We investigate the scaling properties of the Barkhausen effect, recording the noise in several 
soft ferromagnetic materials: polycrystals with different grain sizes and amorphous alloys. We 
measure the Barkhausen avalanche distributions and determine the scaling exponents. In the limit 
of vanishing external field rate, we can group the samples in two distinct classes, characterized by 
exponents r = 1.50 ± 0.05 or r = 1.27 ± 0.03, for the avalanche size distributions. We interpret 
these results in terms of the depinning transition of domain walls and obtain an expression relating 
the cutoff of the distributions to the demagnetizing factor which is in quantitative agreement with 
experiments. 



The Barkhausen noise is an indirect measure of com- 
plex microscopic magnetization processes and is com- 
monly employed as a tool to investigate ferromagnetic 
materials [Q. In recent years, the interest for this phe- 
nomenon has grown considerably due to the connections 
with disordered systems and non-equilibrium critical phe- 
nomena. 

Experiments have shown that the histogram of 
Barkhausen jump (avalanche) sizes follows a power law 
distribution [^-^, suggesting the presence of an under- 
lying critical point |Q,0,D. This hypothesis implies that 
the statistical properties of the noise should be described 
by universal scaling laws, with critical exponents that are 
independent of the material microstructure, at least for 
some class of materials. However, the exact nature of the 
critical behavior is still debated. Several authors |§,||,||j9| 
have analyzed the dynamics of flexible domain walls in 
random media relating the Barkhausen exponents to the 
scaling expected at the depinning transition. Other the- 
oretical explanations involve a critical point tuned by the 
disorder in the the framework of disordered spin models 
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The experiments reported in the literature can hardly 
resolve the theoretical issues since the measured expo- 
nents span a relatively wide range |-|,|ri| and it is difh- 
cult to confirm whether universality holds. In particular, 
there has been no extensive and systematic measurement 
of critical exponents in different materials with homoge- 
neous and controlled experimental conditions and reliable 
statistics. For example, experimental evidence of univer- 
sality has recently been reported for acoustic emission 
avalanches recorded during martensitic transformations, 
by analyzing in detail the behavior of several alloys with 
different compositions [ p^ . 

The precise dependence of the Barkhausen character- 
istic sizes on experimentally tunable parameters is still 
a debated question [p|j6|-|s|JTl[| and a complete agreement 
between theory, simulations and experiments is still lack- 
ing. Understanding this point is crucial in order to link 
the material microstructure to the noise properties, or 
conversely to use the noise to obtain information on the 
structure of the material. In the past, this problem. 



which has important technological relevance, has been 
mainly addressed using phenomenological models p]p^] 
where the material properties (i.e. grain size ||l4(7 in- 
ternal stresses |l^) are accounted for by effective fitting 
parameters. The main limit of this approach is that there 
is no systematic way to implement the program, without 
a precise understanding of the microscopic dynamics. 

In this letter, we report experimental data for a large 
set of materials suggesting the existence of two distinct 
universality classes and show that the results are in quan- 
titative agreement with the theory of domain wall depin- 
ning transition. We perform Barkhausen noise measure- 
ments in six different ferromagnetic materials under sim- 
ilar experimental conditions, averaging the distributions 
over a large number (~ 10^ — 10^) of events, carefully 
testing the effect of the magnetic field rate on the ex- 
ponents. The cutoff of the distributions is tuned by the 
demagnetizing factor and defines two new critical expo- 
nents, which we measure for two materials belonging to 
the different classes. Using scaling arguments, we predict 
values for these exponents that are in good agreement 
with experiments. 

We record the Barkhausen noise using standard induc- 
tive methods, described in details in Refs. |^,||,^. A long 
solenoid provides an homogeneous low frequency trian- 
gular driving field, while a secondary pickup around the 
sample cross section gets the induced flux. The solenoid 
is 60 cm long, with a value of 1450 turns/meter, ensur- 
ing an homogeneous field up to 55 cm long samples with 
peak amplitude of about 150 A/m. The pickup is made of 
50 isolated copper turns, wounded within 1 mm. Such a 
small width is required to avoid spurious effects due to de- 
magnetizing fields. All the measurements are performed 
only in the central part of the hysteresis loop around the 
coercive field, where domain wall motion is the relevant 
magnetization mechanism [|l[. We take special care to 
reduce excess external noise during the measurements of 
avalanches distributions, as the evaluation of critical ex- 
ponents is strongly affected by spurious noise. In this 
respect, the most appropriate cutoff frequency of the low 
pass pre-amplifier filter is chosen in the 3-20 kHz range, 
roughly half of the sampling frequency, as usual in noise 
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measurements. 

The analysis of Barkhausen avalanche distribution is 
performed following the procedure discussed in Ref. d). 
We impose a reference level for Vr for the signal v{t), cho- 
sen above the background noise. The duration T of the 
Barkhausen avalanches is defined as the interval within 
two successive intersections of the signal with the v = Vr 
line. The avalanche size s is calculated as the integral of 
the signal between the same points. We observe that the 
avalanche distributions follow a power law 

p(s) = s-v(.s/so), P(T) = r-"g(r/ro), (i) 

where sq and Tq indicate the position of the cutoff to the 
power law behavior. The critical exponents result to be 
independent of the reference level for a reasonable range 

We employ several different soft magnetic materials, 
both polycrystalline and amorphous: an Fc-Si 7.8 wt.% 
strip (30 cm x 0.5 cm x 60 /xm) produced by plan flow 
casting, annealed several times around 950°C to obtain 
grains of average dimension of 25 /im; two strips of Fe- 
Si 6.5 wt.% (30 cm x 0.5 cm x 45 /itm), one annealed 
for 2h at 1200°C, with grains of 160 /xm, and the other 
annealed for 2h at 1050°C, with grains of 35 /xm [ p^ . 
The amorphous samples have composition of the type 
Fea;Co85_:rBi5 and we employ Fe2iCo64Bi5 as cast (20 
cm X 1 cm X 22 /im), Fe64Co2iBi5 as cast (28 cm x 1 
cm x 23 /im). With these highly magnetostrictivc alloys 
(As ~ 30 - 50 X 10-^) a tensile stress of cr - 100 MPa 
is applied during the measurement. The applied stress is 
found to enhance the signal-noise ratio, reducing biases 
in the distributions, but does not change the exponents 
[ pTf . A partially crystallized Fe64Co2iBi5 (22 cm x 1 cm 
X 23 /im) is also employed, with annealing for 30 min at 
350°C and then for 4h at 300°C under an applied tensile 
stress of 500 MPa. This induces the formation of a-Fe 
crystals of about 50 nm, with a crystal fraction of ~ 5% 

In Fig. la we show the avalanche size distribution, 
obtained for the smallest available magnetic field rates 
(/ = 3-5 mHz). We clearly see that the data can be 
grouped in two universality classes with r = 1.50 ± 0.05 
and T = 1.27 ± 0.03. The first class includes all the Si- 
Fe polycrystals and the partially crystallized amorphous 
alloy, while the amorphous alloys under stress belong to 
the second class. For the materials in the first class, we 
observed a linear decrease of the exponents on the fre- 
quency / of the external magnetic field, in agreement 
with earlier findings Q . The material in the second class 
do not show any noticeable dependence of the exponents 
on the field rate. We note that r ~ 1.3, independent of 
the frequency, was previously measured in Perminvar [Q . 
Next, we measure the distribution of avalanche durations 
(see Fig. lb) and find a = 2.0 ± 0.2 and a = 1.5 ± 0.1 for 
the two classes, with a quite large error bar due to the 
limited range of scaling, and the presence of unavoidable 
excess external noise at low durations. Also in this case. 



a decreases linearly with / for the materials belonging 
to the first class. 

The scaling of the cutoff of Barkhausen avalanche dis- 
tributions has been the object of an intense debate in 
the literature ||6|-|9|,p| . In Ref. [|j the control parameter 
was identified with the demagnetizing factor k. We thus 
measure the Barkhausen avalanche distributions varying 
k, using samples of different aspect ratios. In particu- 
lar, we use the same sample and cut it progressively in 
shorter pieces, recording the noise always in the same re- 
gion, whose size is limited by the pickup coil width. In 
this way only k is varied, while stress, internal disorder 
and system size are kept constant. The demagnetizing 
factor is estimated as A: = ^/^J■c — ^/fJ-i where fic is the 
linear permeability around the coercive field and /i^ is the 
intrinsic permeability (i.e. in an infinite strip) estimated 
using a magnetic yoke jl^ . 

We perform the measurements on the Fe-Si 6.5 wt.% 
1200°C (with lengths spanning from 28 to 10 cm) and 
the Fe2iCo64Bi5 (from 27 to 8 cm) under constant ten- 
sile stress. In Fig. 2 we report the avalanche size dis- 
tribution for different k for Fe2iCo64Bi5 in order to 
show the increase of the cutoff. Data collapse analysis 
yields sq ~ fc~^/°''= with l/tJfc ~ 0.78 (see the inset of 
Fig. 2). Similarly the duration distribution cutoff scales 
as So ^ k~^'' with ~ 0.4. In the case of Fe-Si, the 
analysis is complicated by the frequency dependence of 
the exponents, therefore we fit the cutoff for different 
values of / and extrapolate the results for / ^ 0. The 
results for sq and Tq for both materials are reported in 
Fig. 3 and Table |. 

To interpret the experimental results we use the model 
of domain wall depinning discussed in Ref. |^,|| . A single 
180° domain wall is described by its position /i(r), divid- 
ing two regions of opposite magnetization directed along 
the X axis. The total energy for a given configuration is 
the sum of different contributions due to magnetostatic, 
ferromagnetic and magneto-crystalline interactions, and 
gives rise to the following equation of motion 

^dh(^ = 2noMsH ~k [ (fr'h{f\ t) + -f^V^h{r, t) 



dt 



+ / d^r'j{f- r')ih{f') - h{r)) + T]{f, h), 



(2) 



where F is the viscosity, Ms is the saturation magnetiza- 
tion, H is the applied field, k is the demagnetizing factor, 
7,ii is the surface tension of the wall, J is the kernel due 
to dipolar interactions given by 



J(7--f') 



27r f 



:*/|3 



1 



(3) 



and ?7(r, h) is a Gaussian uncorrelated random field tak- 
ing into account all the possible effects of dislocations, 
residual stress and non-magnetic inclusions. 

The critical behavior of Eq. ^ has been understood us- 
ing renormalization group methods PQ-E2|, which show 
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that at large length scales the critical exponents take 
mean- field values 1^,^. This result is due to the linear 
dependence on the momentum of the interaction kernel 
(Eq. H) in Fourier space |23|. In general, if we consider 
an interface whose interaction kernel in momentum space 
scales as J{q) = Jol?]^, the upper critical dimension is 
given by dc — 2/i and the values of the exponents depend 
on /i (see Table ^. In particular, Eq. |2| yields r = 3/2 
and a = 2 (i.e. = 1), or t ~ 1.27 and a ~ 1.5 if 
dipolar interactions are neglected (i.e. n = 2) Nu- 
merical simulations confirm the linear dependence of the 
exponents on the driving frequency / for /i = 1, and no 
dependence for = 2 ||l^ . 

The experimental results are in perfect agreement 
with the values predicted using Eq. || and suggest that 
the dipolar interactions are stronger than surface ten- 
sion effects in polycrystals, or whenever small grains 
are present, while in amorphous alloys under stress 
the surface tension is much stronger. Magnetostric- 
tion can be one of the sources of this effect, since the 
domain wall surface tension increases with stress a as 
7u, ^ ^ Kq + 3/2AsCr, where Kq is the zero applied 
stress anisotropy and is the saturation magnetostric- 
tion [p]p7|]. Micromagnetic calculations for these partic- 
ular materials are needed to have a final confirmation of 
this effect. 

We derive the dependence of the cutoff on k from Eq. |^ 
noting that the demagnetizing field acts as a restoring 
force for the interface motion and is responsible for the 
cutoff in the avalanche sizes. The interface can not jump 
over distances larger than ^, defined as the length for 
which the interaction term Jo|g|'' is overcome by the 
restoring force (i.e. JohS,~'^ ~ kS^'^h) which implies 

i-{klJo)-'"' i/fe = l/(M + d). (4) 

The avalanche size and duration distributions cutoff 
can be obtained using the scaling relations reported in 
Ref. §: 

sor^D{k/Jo)-^'''\ + (5) 

and similarly 

To^D{k/.h)-^'^, Afc = z/.fc, (6) 

where D = ^ {rf) denotes the typical fluctuation of the 
disorder. The dynamic exponent z and the interface 
roughness exponent C, define the spatio-temporal scaling 
of the domain wall width ^ {K^) - {h}^ = ^'^ -F(i/D- 

Inserting in Eqs. (|^^ the renormalization group re- 
sults C = (2Ai - rf)/3 and z = ^ - (4^ - 2d)/9, we 
obtain l/at = 2/3 and = (/^ - (4// - 2d)/9)/(^ + d). 
We have performed extensive numerical simulations in 
d — 2 using the model described in Ref. [|| in order to 
test these results (see Table Q). Our results also agree 
with earlier numerical simulations in d = 1 for 1 < /i < 2 
where sq ~ independent of /i |24] and with the re- 
sult reported in Ref. [ni| (i.e. sq ^ Lr^'^) obtained using 
k = l/L\ 



The scaling of the cutoff is different for a local demag- 
netizing field —kh{x,t) [Q, which yields i/fc = l//i 
However, a non-local term is more appropriate to de- 
scribe the demagnetizing field, due to the long-range of 
magnetostatic interactions Q, as it is confirmed by the 
agreement between experiments and theory. 

In conclusions, our experiments suggest that the 
Barkhauscn effect can be described by universal scaling 
functions and that materials can be classified in different 
universality classes. The theory of interface depinning 
can be used to obtain a quantitative explanation of the 
experiments, providing a natural framework to under- 
stand the properties of soft magnetic materials. 
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1/afc 




RG 


3/2 
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2/3 


1/3 


Sim 


1.50±0.02 


2.00 ±0.05 


0.65±0.05 


0.3±0.1 


Exp 


1.50±0.05 


2.0±0.2 


0.57±0.10 


0.30±0.10 


^1 = 2 


r 


a 


1/afc 


Afc 


RG 


5/4 


10/7 


2/3 


7/18 


Sim 


1.27±0.02 


1.50 ±0.05 


0.72±0.03 


0.39±0.05 


Exp 


1.27±0.03 


1.5±0.1 


0.79±0.10 


0.46±0.10 



TABLE I. The critical exponents measured in experiments 



compared with RG results 



23| valid to order € = 2/1 — d 
(see Ref. |6| for the scaling relations) and simulations for the 
two universality classes. System sizes range up to L = 150 
for /J. = 1 and to L = 400 for jj, — 2. 
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FIG. 1. (a) Distributions of Barkhausen jump sizes mea- 
sured in different materials for the lowest available driving 
frequency. The solid line has a slope r = 1.5 while for the 
dashed one r = 1.27, corresponding to the two universality 
classes; (b) Similar plot for duration distributions. The solid 
line has a slope a = 2 while for the dashed one a = 1.5. 
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FIG. 2. The Barkhausen jump size distribution for differ- 
ent values of the demagnetizing factor k in the Fe2iCo64Bi5 
amorphous alloy under tensile stress The data collapse re- 
ported in the inset is done using r = 1.27 and 1/ak ~ 0.78. 
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FIG. 3. The cutoff of the Barkhausen size (a) and duration 
(b) distributions as a function of k in Fe-Si 1200° C (empty 
symbols) and Fe64Co2iBi5 under constant tensile stress (filled 
symbols) . 
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